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Introduction. A three-dimensional dry friction model in the interaction of a rectangular body and a horizontal rough 
surface is considered. It is assumed that there is no separation of the body from the horizontal surface. The body motion 
occurs under the conditions of combined dynamics when, in addition to the longitudinal movement, the body partici- 
pates in twisting. 

Materials and Methods. Linear fractional Pade approximations are proposed, which replaced the cumbersome analytical 
expressions that most accurately describe the motion of bodies on rough surfaces. New mathematical models describing 
sliding and twisting of bodies with a rectangular base are proposed. 

Results. Analytical expressions of the principal vector and moment of friction for rectangular contact areas are devel- 
oped and scientifically established. A friction model that takes into account the relationship between sliding and twist- 
ing speeds, which provides finding solutions for Pade dependences, is developed. After numerical solution to the equa- 
tions of motion, the dependences of the sliding speed and angular velocity on time were obtained and constructed. 
Graphs of the dependences of the friction forces and their moment on two parameters (angular velocity and slip veloci- 
ty) were constructed, which enabled to compare the integral and normalized models of friction. The comparison results 
showed good agreement of the integral model and the model based on Pade approximations. 

Discussion and Conclusions. The results obtained provide considering the dynamic coupling of components, which 
determines the force interaction of a rectangular body and a horizontal surface. These results can be used in mobile ro- 
botics. The analyzed motion of the body occurs through the motion control of a material point inside the body. Such 
mobile robots can be used when solving a wide class of problems: when creating autonomous robots for the exploration 
of outer space and planets; in the diagnosis and treatment in case of passing through complex structures of veins and 
arteries; in research under water, in places of large differential temperature; in underground operations. 
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Introduction. The study of the movement of a rectangular body is a challenge in the mobile robotics [1]. This 
movement is due to the control of the material point inside the body. Such mobile robots can be used to solve a wide 
range of tasks. For example, when creating autonomous robots for the exploration of outer space and planets; for 
medical purposes, in diagnosis and treatment, for example, in case of passing through complex structures of veins and 
arteries; as well as for underground work and research under difficult conditions, for example, under water and in places 
of large differential temperature [1, 2]. 

Thus, more and more challenges are being set for robotics, which require theoretical research, including 
studying models of friction between the body and the surface under the conditions of combined dynamics [3, 4]. Since 
the movement of the mobile robot occurs in different directions, it is required to consider the longitudinal movement 
and rotation. Thus, in the structure of the friction model, it is required to provide the relationship between the sliding 


and twisting speeds [5]. An important development in the description of this relationship was made in [6]. Its author 
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managed to solve the equations for the principal moment and the vector of friction forces where a rectangle was 
considered as a contact area. Such analytical expressions enable to most accurately describe the motion of bodies on 
rough surfaces, but they are cumbersome and complex since they contain integral expressions. Hence, the authors of [7] 
constructed linear fractional Pade approximations, which made it possible to find solutions for the resulting 
dependences. 

Pade approximation can be used to explain the effects of combined dry friction for linear and angular 
velocities. On the basis of Pade approximations, it became possible to create new models of friction [8, 9], which later 
began to be classified for better interpretation [10]. The classification occurs depending on the number of parameters. 
Thus, in [11], the authors introduced the notions of dimension and order of the dry friction model depending on the 
order of the used Pade approximations. 

The model of sliding and twisting friction, which is proposed in the paper [12], provides considering the 
dynamic connection of the components that determine the force interaction of a rectangular body and a horizontal 


surface [13]. 


Problem Statement. We consider a solid body of mass m, , which is a rectangular body with uniform faces of 
length a, width 5 and height 2h. A fixed coordinate system Oxyz , associated with the body (Fig. 1) is introduced. 
Point O is located on the horizontal plane. The system. O,x,y,z, starts at a point O,, that corresponds to the geometric 
center of the body. The axis Oz, is parallel to the axis Oz. The axis O,x, 1s parallel to the long edge of the body. We 


introduce the unit vectors e.,e, of the axes O,y, and O,x,, respectively. 
were 








Fig. 1. Coordinate systems 


Consider the continuous motion of the body on the surface (Fig. 2), which consists of translational movement 
and rotation about the axis O,z,. Three coordinates determine the position of the body. The coordinates x,, y, and A 


set the origin of the coordinate system O,x,y,z, in the coordinates Oxyz . The rotation of the body relative to its initial 
position on the axis O,x, is specified by the angle @. This paper considers the case when the center of mass of the body 


G and the center of mass of the system O, coincide (Fig. 2) [14]. 
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xy 





mg 


Fig. 2. Movement of the system body 
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Materials and Methods. The contact area is a rectangle with sides a and b, in which the normal voltage 
depends on the distance from the point P to the faces of the rectangle (Fig. 3). 





Fig. 3. Velocities of the points O, and P 


Consider an infinitesimal area dS at an arbitrary point M on the contact surface. We introduce the angle 
between the relative sliding velocity and the axis O,x,. Let us draw the radius vector 7,,, from the point P to the point 


M. The velocity vector of the point M is denoted uy, and to find it, we use the Euler formula describing the velocity 
distribution in a perfectly rigid body: 
Uy = Up + W-Typ. 
The sliding speed at the point M is decomposed into two components along the axes 0,x, and O,y,: 
Umx = Vx — YO; 
Umy = Vy + XW. 
Using Coulomb's law, we find a small increment of the friction force directed against the relative velocity at 
the point M [15]: 
dF = —fo(x, yids, 


where f — coefficient of friction; o(x, y) — contact stress distribution function depending on the x and y coordinates; 


dS = dxdy — small area increment [15]. 


We rewrite the differential of the friction force and the moment of this force in projections on the axes under 





consideration: 
dE, = —fo(x, y)—™ dxdy; 
luv 
dF, = —fo(x,y) -“Y dxdy: 
lum 
i jk 
dM,=|*xX yy O|=xdF, — ydF,. 
dk, dF, 0 


As a special case, we consider a uniform distribution of stress in the absence of internal masses in the body, 
then these stresses will be equal to: o = a but then we will continue the record in general form: o(x, y). 


Having integrated the expressions for the friction forces, we obtain: 








/2 pb/2 = 
F, = —f ‘en | 5/3 o(x,y) ar dxdy; (1) 
a/2 -b/2 Vytxw 
Fy a = ine Spyz o(x,y) i dxdy. (2) 


Relative slip module |v,,| is calculated from the formula: 


luyl = v3, + Viry = v2 gO (ear y-) 2w(v,x — Uxy). (3) 
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Imagine the relative positions of the vectors of variable sliding speed uv and the components of the friction 
force: F — the component opposite to the sliding speed v; F,; — the component perpendicular to the instantaneous slip 


velocity. At the same time, imagine the coordination of this system with respect to the axes O,x, and O,y, (Fig. 4). 





Fig. 4. Components of the friction force and velocity 


We will make a transition from the projections of the sliding speed: 


Vy = UCOS®, 
is = vsing, (4) 


to the speed module and the sliding angle: 


F, = F,cos@ + F,sing, 


F, = F,(-sing) + Fsing’ (5) 


We will integrate the moment of the friction force on the contact area: 
a/2 (gee u(xsing—ycos@)+w(x2+y7) 


M, = —f Jays —p/2 a(x, y) xdy. 


v2+02(x2+y2)+2wvu(xsing—ycos@) 


Let us substitute the expressions (1)—(3) into the system (5), and also rewrite the expression for the moment of 
force. As a result, we obtain a three-dimensional model of friction sliding and twisting: 
a/2 


b/2 v(cos?@+sin? ~)—w(ycosptxsing) 


Py cae’ d lea ye o(x, y) v24+w2(x2+y2)+2wv(xsin~—ycos@) (6) 
= a/2 b/2 w(ysin@+xcosq@) 

P= poe? Son72 o(x, y) v24+w2(x2+y2)4+2wv(xsinp—ycos@) 7) 
_ _ ¢ fa/2 pb/2 u(xsing—ycos@)+w(x2+y7) 

M, = Pye Son72 o(x, y) v24+w2(x2+y2)+2wv(xsing—ycos@) aay co 


In order not to solve cumbersome integrals, we use the replacement of the corresponding Pade expansions [16, 
17]. Thus, based on the Pade theory [18], these expressions can be formulated as the ratio of two functions of several 
variables in the entire domain of definition, provided that the functions must have the same order [7]. To define these 


functions, it is required to determine the behavior of integral expressions (6)—(8) under the following conditions: 











OF | _ _ f a/2 pb/2 x2 +y?(x*+y7)+(ycosptxsing) (xsinp—ycos@) =: oy. 
dv \v=0 7 w ?—-a/2 ee a(x, y) (x2+4y2)2 dxdy ~ Ww Io; 
aM, _ _ f pa/2 b/2 _ 2y?xsing—2y* cos@) a ae 
dv lu=0 wt -a/2 pe o(x, y) (x2 4+y2)2 dxdy 133 
a/2 ¢b/2 x2 —y? 

Mo e900 a =f Joss eee o(x,y) " Were iad = = fle. 

a/2 b/2 ysin~+xcos@ 

Figoo = —F Saj2 Sp). TY) Tlgye  OXAY = flo; 
oh set ser 2 rr ae 
OV jo=0 IW|y=0 IW|w=0 Nese OV Jw=0 IW |y=0 

_ _ OF, _ OM, _ OM, _ _ 

—_ Pies ~ da l=0 ~ Oy io=6 Ba lv=0 _ M7 1y-00 = 0. 
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. OF aM 
Values of the expressions —= = 


are not involved in finding the subsequent Pade approximants, 
dV |v=0 dW |w=0 





therefore, their writing is omitted due to their cumbersomeness. The identical equality to zero is realized under the 
condition that the voltage o is symmetric about the center of the rectangular contact spot, 1.e., the point P. 

An accurate three-dimensional integral model [13] (6)-(8) provides a logical description of dry friction 
phenomena, but for solving problems of dynamics, such a model is difficult to accept due to the need to calculate 
impressive integrals [10]. To avoid this procedure, we use [6] to replace the exact integral system with the 
corresponding expressions using Pade approximations in the whole range of variables. Linear fractional Pade 


expansions give a three-dimensional model of first-order sliding and twisting friction [19]: 











v+b1w 
Fi = Fo ao (9) 
_ W+b2v_ 
M, = Mo (10) 
_ W+b3u 
Td aera (11) 


To determine the Pade coefficients, it is required to study the properties of this model at the boundary points 
by analogy with integral expressions. To do this, we differentiate the parameters F), F,, M, and thus satisfy the 


corresponding integral expressions: 


(Fai —— 


fo 
vt T, WO 
M, = —flg 
I OO 
Fy = -fle—— 
{ w+ Fv 
9 
The system of equations of motion has the form: 
aWo, 
J dt Z 
dV x 
(mp +m) _ = F, + (mp + My, )VyWo,; (12) 
dvy 
(Mm, + m,) “ae _ ky —- (mM, at M)Vx.Wo,- 
We express the time derivatives of the sliding speed and the angular speed using the formulas (3)-(5): 
du 1 ; dv, 5 dv, 
— = ———— | 20, —— + 20, — |; 
dt 2./v2 +02 ~ dt Y dt 
d 1 v sin v cos 
a (Bain ens 
dt m)j+m, vu? v2 


We rewrite these equations using the formulas (3)-(5) for Wo, = w/a and add the first equation from the 
system (12): 
J ad = M,a; 
dt — 
(Mp + my) < = Fy; (13) 
(m,) + m,)Uuq@ = F,. 


Research Results. Next, we calculate the integral expressions of the parameters Jp, 1,, Iz, Ig, Ig c using the 


Wolfram Mathematica software package for the following values: 


f=1,a= 05m; b= 0.2m; my = 1kg, o= a a ee 
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and substitute in the system of equations (13). Based on numerical expressions, we build graphs of integral and 
normalized functions depending on the parameter = ¥/,, . Fig. 5 shows graphs of the functions of the integral friction 


els (11)-(13), as well as models based on Pade approximations (14)-(16). 
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a) k b )) k 
Fig. 5. Graphs of integral (solid lines) and normalized (dotted lines) functions of the tangent (a), normal (b) components of friction 
force and friction moment (c) 


Based on the graphs of the functions (Fig. 5), we can talk about good matching of the considered models. 


Next, we obtain graphs of the dependences of the characteristic parameters on time (Fig. 6). 
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Fig. 6. Dependences of sliding velocity v (a) and angular velocity w (b) on time ¢ 


Discussion and Conclusions. The movement of the mobile robot, starting from the contact of its body and the 
reference plane, under the conditions of combined dynamics, when there is sliding and twisting, is described. Analytical 
integral expressions are obtained for the tangent and normal components of the friction force [19] and the moment of 
friction applied to a rectangular contact area. The corresponding Pade approximations are determined for the obtained 
expressions. The integral and normalized models are compared through plotting the dependences of the friction forces 
and the moment of friction on the angular velocity and the slip velocity. The comparison results showed good matching 
of the integral model and the model based on the Pade decompositions. The graphs correspond to the logical behavior 
when a rectangular body moves, since the sliding speed and angular velocity increase according to the specified 
parameters. Consequently, the combined friction model implemented using Pade approximations can be applied to 
solve problems related to mobile robots with a rectangular base. 
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